PROCEEDINGS 

°F SCIENCE 



00 

o 
o 

(N 
D 
OO 



(N 
> 

in 

00 
OO 
(N 

G\ 

o 

OO 

o 



Nonperturbative infrared fixed point in sextet QCD 



Benjamin Svetitsky* and Yigal Shamir 

Raymond and Beverly Sackler School of Physics and Astronomy, Tel Aviv University, 69978 Tel 
Aviv, Israel 



E-mail: bqs@julian.tau.ac.il shamir@post.tau.ac.il 



Thiomas DeGrand 

Department of Physics, University of Colorado, Boulder, CO 80309, USA 
E-mail: [degrand@pizero . Colorado, edu 



The SU(3) gauge theory with fermions in the sextet representation is one of several theories of 
interest for technicolor models. We have carried out a Schrodinger functional (SF) calculation 
for the lattice theory with two flavors of Wilson fermions. We find that the discrete beta function 
changes sign when the SF renormalized coupling is in the neighborhood of g-^ — 2.0, showing a 
breakdown of the perturbative picture even though the coupling is weak. The most straightforward 
interpretation is an infrared-stable fixed point. 



The XXVI International Symposium on Lattice Field Theory 
July 14-19, 2008 
Williamsburg, Virginia, USA 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Nonperturbative infrared fixed point in sextet QCD 



Benjamin Svetitsky 



1. Motivation — Beyond the Standard Model 

Most interesting theories that go beyond the Standard Model require nonperturbative informa- 
tion in order to demonstrate their relevance to low-energy physics Among these are strong- 
coupling realizations of the Weinberg-Salam theory; technicolor, in which the Higgs multiplet 
arises as bound states of a higher-energy theory; and theories based on extra dimensions, in which 
Kaluza-Klein phenomena may depend on strong-coupling physics. Moreover, various theories of 
unification start with large gauge groups that reduce to low-energy theories via tumbling and vac- 
uum alignment, phenomena that are inherently nonperturbative. Supersymmetry can only apply to 
the low-energy world if it is broken by some nonperturbative mechanism. 

As a first step in what we hope will develop into a broad attack on gauge theories beyond QCD, 
we have chosen to study the SU(3) gauge theory with two flavors of Wilson fermions in the sextet 
representation [|2|]. The two-loop beta function of this theory crosses zero Igll] at^2^ 10.4, which 
is a strong coupling; a ladder calculation indicates that the quarks condense and chiral symmetry 
is spontaneously broken before this coupling is reached. If indeed chiral symmetry is broken in 
this theory, it becomes a candidate for a theory of walking technicolor [^], assuming that a lattice 
calculation can confirm the slow evolution of the coupling. The survival of the fixed point, on the 
other hand, would put the theory in the conformal window and disqualify it for technicolor. 

A related issue is the possibility of scale separation, where the confinement scale of the theory 
is at a lower energy than the scale of the chiral condensate |^]. This is what initially attracted us to 
studying quarks in higher representations than the fundamental. 



2. Perturbative renormalization group 

Let me begin by reviewing the possibilities raised by the two-loop beta function, as described 
by Banks and Zaks [J]. The perturbative expansion is 

where b\>0 and b2 < 0. If we truncate at these two terms, the formula give an IR-attractive fixed 
point (IRFP) at g = g^,, as shown in Fig. [T[ 

If is a weak coupling, then the two-loop calculation may be reliable. Flow into the IRFP 
implies that the massless theory possesses conformal dynamics at large distances. This means that 
there is no confinement, no chiral condensate, and indeed no particles whatsoever!^ 

If g^ is a strong coupling, on the other hand, the chiral condensate will form before the IRFP is 
reached, so the quarks will become massive and decouple from the IR dynamics. The beta function 
past that point, returning to that of the pure gauge theory, will remain negative all the way to ^ = oo 
and there will be no actual zero. The sketch for this case in Fig. |] shows the marginal possibility 
that the beta function hovers near zero as this decoupling takes place, which is what underlies the 
hypothesis of walking. 



'We refer here to Casher's argument [Q] that massless quarks cannot form bound states in a vector-coupled gauge 
theory unless chiral symmetry is spontaneously broken. All color-singlet Green functions will possess only cuts starting 
at = 0. 
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Figure 1: Sketch of the beta function if the two-loop zero occurs at weak coupling (left) or at strong 
coupling (right). The green arrows on the axis are the directions of the IR flow. The shape of j3 [g) on the 
right indicates walking. 

The two-loop beta function makes a prima facie case for the existence of an IRFP and thus 
places the theory in what is called the conformal window [||]. As reviewed in [||], however, a calcu- 
lation with the Bethe-Salpeter equation (the "ladder" approximation) points to ;tSB; a conjecture 
based on supersymmetry on the other hand, puts the theory back in the conformal window. A 
lattice determination of the beta function can resolve the matter. 

3. The Schrodinger Functional method 

The Schrodinger functional (SF) [|lO[ |l|] 

is a well-known method for calculating the beta 
function of the theory via imposing a background field. Following the method of SF calculations 
for QCD, we employ Wilson fermions rather than staggered so that boundary values U[-, (the back- 
ground field) can be set on single time slices at f = and ? = L on an lattice; Wilson fermions 
also give us better control over the number of fermion flavors. We add a clover term to remove 0{a) 
discretization errors, and we fix the clover coefficient self-consistently via tadpole improvement. 

Wilson fermions, of course, break chiral symmetry explicitly and one must fix K" = jq. to have 
a massless theory in the continuum limit. We define the quark mass by the axial Ward identity, 

_1 d,{Al{t)(f\t' = 0,p = 0)) 

where if is an operator on the boundary at ? = while A4 and P are the axial and pseudoscalar 
densities, measured at zero spatial momentum at the center of the lattice. Tuning to niq = fixes 

As is usual in SF calculations, we give the fermion fields a spatial twist. The boundary con- 
ditions then serve as an efficient IR cutoff, even stabilizing the fermion inversions at = jq. and 
allowing us to study the massless theory directly. 

In the continuum limit, the background field depends only on the size L of the system, so the 
method gives the running coupling at the IR scale, g^{L). More precisely, one is to calculate the 
potential F = — logZ and compare it to the classical Yang -Mills action 5^' of the background field 
configuration, giving g^{L) via 




(3.2) 
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In a lattice Monte Carlo calculation, however, one cannot calculate T directly. The trick is to let the 
boundary values U' depend on a parameter T]. Then by differentiating Eq. (3^) we relate g^{L) to 
the expectation value of an operator, viz.. 



5r 



dS. 



gauge 



tr 



1 d{DlDF) 1 



K 



K 



37.7. 



(3.3) 



To summarize: In order to extract the scaling of the running coupling, we 

1. fix the lattice size L and the couplings p and K = JQ-(j8); 

2. calculate K/g^{L) via the expectation value (33), and also 



3. calculate K/g^(2L) on a lattice twice as large. The two lattices have the same bare parameters 
(j8, k:) and hence the same UV cutoff a. We thus 



4. obtain the discrete beta function (DBF), defined as the difference 

. , K K 
B{u,2) 



(3.4) 



g^i2L) g^{L) ' 
which is a function ofu = K/g^{L). 
We show the result of this procedure in Fig. ^(left). The DBF B{u,2), obtained from lattices 
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Figure 2: Left: Discrete beta function B{u,2) obtained by comparing lattices of size 4"* and Right: 
B(m,4/3) obtained from lattices of size 6'* and 8^. The dashed curves are the two-loop predictions. 

of size 4^ and 8^, evidently crosses zero at g^ ~ 2.0. This is a much weaker coupling than ~ 10 
as found in two-loop perturbation theory. (For comparison, the two-loop DBF is plotted as a dashed 
curve.) On the face of it, this result demonstrates that the massless theory possesses an IRFP and 
hence that the IR theory is conformal. 

If a lattice of size 4^ seems small, we can compare instead lattices of size 6^ and 8^, which 
yield the DBF S(m,4/3) for the smaller scale factor 4/3 (Fig. ^ right). As one might expect, the 
result is generally closer to zero than B{u,2), and hence the error bars are relatively larger. The 
crossing of zero is evident nonetheless. We note that the data points in each plot are statistically 
independent, but the two plots are linked since the same 8"^ data are used in each. 
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4. Caveat cursor 

The DBF is, in principle, a continuum quantity that relates couplings at different IR scales L. 
A lattice calculation of the DBF introduces an implicit dependence on the lattice spacing. It goes 
without saying that our results contain lattice artifacts, which we cannot estimate since we have 
worked so far at a single lattice spacing for each value of g and for each rescaling factor (2 and 
4/3). 

Even when a satisfactory continuum limit is reached, however, one must ask whether the 
theory is really described by a single running coupling. After all, in the IR regime any theory will 
need more than one term in its Lagrangian to describe its spectrum and interactions, and these terms 
will be generated by RG transformations. Continuum perturbation theory automatically limits the 
effective Lagrangian to renormalizable couplings, and thus one always speaks of a single beta 
function for QCD or other (massless) gauge theories. A lattice theory, on the other hand, can 
quickly generate many terms in the effective Lagrangian unless a physically reasonable truncation 
is used in the RG transformation. 

One necessary condition that the lattice theory be well described by a single coupling at scale 
L is to verify that L is too small for confinement to have set in. Figure |3| is the phase diagram of the 
lattice theory. It is clear that on our lattice the IRFP is found in the weak-coupling phase, meaning 
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K 

0.15 
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Figure 3: Phase diagram in the (/3,fc) plane. The solid curve is 'Q ()3), where the AWI quark mass niq 
vanishes. The star marks the couplings corresponding to the IRFP of the SF effective coupling. The other 
curve denotes the finite-temperature/finite-volume phase transition for L — %. Presumably the finite-volume 



curve intersects the curve, but we have not gone there. See T. DeGrand's talk 1 12 1 for more information 
on the phase diagram. 

that the IR scale L = 8a is well within the confinement radius (if any). 

If we examine this more closely, the simplicity of a gauge theory at a given scale may be judged 
by the behavior the qq potential. If the potential is almost Coulomb, meaning V{r) = g^{r)/r with 
a coupling g that varies only slowly with r, then one may proceed as if g is the only coupling. 
Again, we expect that this breaks down at large distances in QCD, where V{r) first becomes linear 
in r as the confining flux tube forms and then decays exponentially as it breaks. It may be seen 
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in Fig. ^ that when we measure V{r) we find it to be consistent with a Coulomb potential in the 
neighborhood of the couplings corresponding to the IRFP. 
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Figure 4: Examples of the qq potential from our data. On the left, the string tension is large; the Sommer 
parameter ro/a is measurable and we can obtain a good fit to V{r). On the right, the string tension shrinks 
away as the temporal extent of the Wilson loop grows, and we cannot perform a reliable fit to V{r). The plot 
on the right represents the situation near the fixed point in Fig. |^. The lattice sizes are 8^ x 12 (left) and 12^ 
(right). 

Confirmation of the IRFP will thus come from (1) checking the DBF with more and larger 
volumes, (2) further understanding of the phase diagram in the (j8,K") plane, and eventually (3) 
understanding of the continuum limit. All this is in progress. Then the challenge will be to measure 
properties of the conformal theory at the fixed point. This will entail at least the calculation of 
operator exponents, which will govern how the fixed point is approached from the theories in its 
basin of attraction and from nearby massive theories. 

More work on this model has been presented by Daniel Nogradi at this conference |[l3|]. This 
work was supported in part by the Israel Science Foundation under grant no. 173/05 and by the US 
Department of Energy. Our computer code is based on version 7 of the publicly available code of 



the MILC collaboration |14]. 
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